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ABSTRACT 

Time-distance helioseismology has shown that f-mode travel times contain 
information about horizontal flows in the Sun. The purpose of this study is to 
provide a simple interpretation of these travel times. We study the interaction of 
surface-gravity waves with horizontal flows in an incompressible, plane-parallel 
solar atmosphere. We show that for uniform flows less than roughly 250 ms _1 , 
the travel-time shifts are linear in the flow amplitude. For stronger flows, per- 
turbation theory up to third order is needed to model waveforms. The case of 
small-amplitude spatially-varying flows is treated using the first-order Born ap- 
proximation. We derive two-dimensional Frechet kernels that give the sensitivity 
of travel-time shifts to local flows. We show that the effect of flows on travel 
times depends on wave damping and on the direction from which the observa- 
tions are made. The main physical effect is the advection of the waves by the 
flow rather than the advection of wave sources or the effect of flows on wave 
damping. We compare the two-dimensional sensitivity kernels with simplified 
three-dimensional kernels that only account for wave advection and assume a 
vertical line of sight. We find that the three-dimensional f-mode kernels approx- 
imately separate in the horizontal and vertical coordinates, with the horizontal 
variations given by the simplified two-dimensional kernels. This consistency be- 
tween quite different models gives us confidence in the usefulness of these kernels 
for interpreting quiet-Sun observations. 



Subject headings: convection — scattering — waves — Sun: helioseismology - 
Sun: oscillations 



1. Introduction 

One of the tasks of time-distance helioseismology (Duvall et al. 1993) is to map the flow 
field in the solar convection zone at the highest spatial and temporal resolutions possible. 
Using a cross-correlation technique, it is possible to estimate the time taken by waves to 
travel from one point on the solar surface to another. These travel times contain information 
about local flows (e.g. Duvall & Gizon 2000; Gizon et al. 2000). Particularly interesting are 
flows associated with supergranulation, active regions and sunspots, and connected with the 
evolving magnetic cycle of the Sun. The devil is in the details of how the flow information 
is encoded in and extracted from the travel-time measurements. 

We study the forward problem in the case of solar surface-gravity waves, also called f 
modes, which propagate horizontally. For frequencies near 3 mHz the f-mode kinetic energy 
peaks about 1 Mm below the solar surface. The interaction of f modes with flows has been 
studied by Murawski (2000), Duvall & Gizon (2000), Woodard (2002), Zhao et al. (2007), 
Birch & Gizon (2007), and Birch et al. (2007). Small-scale turbulent convection (granulation) 
excites and damps solar oscillations (e.g., Murawski 2000). Birch & Gizon (2007) and Birch 
et al. (2007) used the Born approximation (e.g. Gizon & Birch 2002) to compute the effects 
of the advection of waves by small-amplitude steady flows on time-distance and ring-diagram 
measurements. 

Here we are interested in studying in more detail the sensitivity of f-mode travel times, 
as measured by time-distance helioseismology, to horizontal steady flows. We wish to answer 
the following questions, none of which have been studied before. How do travel-time shifts 
depend on travel distance? In particular, what is the role of wave damping? Do travel-time 
shifts depend linearly on the flow amplitude for typical solar flows? If so, when does the 
linearity break down? Which physical effects contribute to the travel-time perturbations? Is 
advection of the waves by the flow the only important effect? Or is the Doppler shift of wave 
sources (granulation) by larger-scale flows also relevant? Is it important to model damping 
in a frame moving with the flow? In the limit of weak flows, we can obtain kernel functions 
that give the linear sensitivity of travel-time perturbations to spatially-varying horizontal 
flows. How do sensitivity kernels change with the bandwidth and the central frequency of 
the f-mode wave packets? Waves at the solar surface are typically observed in line-of-sight 
velocity using Doppler measurements. What effect does the projection of the wave velocity 
onto the line of sight have on the kernels? Can we find an accurate approximation to the 
kernels to speed up demanding calculations? Is it possible to recover the ray approximation 
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of Kosovichev & Duvall (1997) from the high-frequency limit of the cross-path integral of 
the sensitivity kernels (see Tong et al. 1998)? 

In order to answer these questions, we choose to work with a simplified model of surface- 
gravity waves. Following Gizon & Birch (2002), we consider waves at the free surface of 
a constant-density half space under constant gravity. We choose initial conditions such 
that the wave field remains irrotational. These assumptions are useful as they simplify the 
presentation; they are also sensible as will be shown at the end of the paper through a 
comparison with a solar-like three-dimensional kernel calculation with otherwise simplified 
physics. 

The layout of the rest of the paper is as follows. In § 2 we introduce our simple surface- 
gravity wave model. In § 3 we study the effect of a uniform flow on wave travel times, finding 
an exact solution and one using perturbation theory. In § 4 we use the first-order Born 
approximation to obtain the wave field in the case of a small-amplitude spatially-varying 
flow, and then derive the sensitivity kernels. § 5 is a comparison with a three-dimensional 
kernel calculation. The implications of our findings are discussed in § 6. 



2. A simplified model for f modes 

Following Gizon & Birch (2002), we consider a constant-density, incompressible half 
space z < 0, bounded by a free surface at z — 0. Throughout the paper we denote the 
position vector by r = (x,z), whose horizontal component is x = (x,y) and z is height. 
Gravitational acceleration, — gz, is constant with g = 274 ms~ 2 . We use a plane-parallel 
geometry as an approximation to a local patch on the Sun: this is valid for wavelengths much 
less than the solar radius. We denote by £ the unit vector pointing toward the observer, 
which depends on the position on the solar disk but is assumed to be independent of position 
over the small solar region of interest. The line-of-sight vector t = (£h,£ z ) decomposes into 
a horizontal component £^ and a vertical component £ z . 

The velocity field has two components: the propagating part, v(r, t), and the underlying 
steady horizontal flow, u(r). As mentioned in the introduction and in order to obtain a 
simple problem, we choose both u and v to be irrotational at all times. We introduce the 
wave velocity potential, 6, such that 

v = VQ. (1) 

The main advantage of this setting is that the problem of surface-gravity wave propagation 
through flows reduces to a two-dimensional one. 

As is appropriate for solar waves, we consider small-amplitude (linear) waves. The 
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linearized hydrodynamic equations we wish to solve are as follows. In the bulk z < 0, the 
momentum and continuity equations are 

-iuQ + u-ve = -p/p-r[e], (2) 
v 2 e = 0, (3) 

where p is the pressure perturbation, p is the density, and T is a phenomenological frequency- 
dependent damping operator (see Gizon & Birch 2002). The above equations and the ones 
that follow apply to quantities that were Fourier-transformed in time, with u the angular 
frequency (we use the Fourier conventions of Gizon & Birch 2002). At the free surface surface 
z — 0, the linearized dynamical and kinematic boundary conditions (e.g., Whitham 1974) 
are 

p/p-grj = U/p z = 0, (4) 
— iurj + u ■ Vr? = d z Q z — 0, (5) 

where 77 is the elevation of the fluid's surface, and n is a stochastic pressure source on the 
surface that generates the waves (Gizon & Birch 2002). 

In many cases, observations of solar oscillations are Doppler velocity observations (e.g. 
Scherrer et al. 1995). Let us introduce the line-of-sight component of the wave velocity at 
the surface and at time t, (f>(x,t), given by 

(f)(x,t)=£-v(x,z = 0,t). (6) 

The actual measurements involve convolution of with the point-spread function of the 
telescope. In addition, the f-mode data analysis includes filtering the observations to select 
the f-mode ridge in Fourier space (Duvall & Gizon 2000). We denote the resulting data cube 
by ip(x,t). The combined effects of the point-spread function and the data filtering can be 
described by the convolution of a function F(x,t) with <fr(x,t): 

il>{x,t) = {F*<j>)(x,t). (7) 

In Fourier space, the filtered data ip(k, u>), where k is the horizontal wave vector, is given by 
a multiplication of the filter F(k ) uo) by 0(fc,cj), 

i/>(k,w) = F(k,u)(f)(k,u). (8) 

As in Gizon & Birch (2002), a mode-mass correction is also included in F to simulate the 
response of a solar-like stratified atmosphere. Note that we use the same symbol for a 
function and its Fourier transform: the arguments are used to distinguish the two. 
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We are first interested in the filtered oscillation signal ip that would be observed in the 
absence of any flows (zero-order terms are denoted with a superscript 0). Eliminating the 
pressure and surface elevation terms, the surface boundary conditions at z = reduce to 

(d z - k)G° = — n = S°, (9) 

pg 

where k{uj) = {u + i^)uj / g is the complex wavenumber at resonance. The zero-order damping 
operator is such that r°[6°] = 76 with the damping frequency ^(uj)/2-r = \uo/uo*\ iA x 
100 /zHz, and u;*/2ir — 3 mHz is the central frequency (Gizon & Birch 2002). 

The solution of equation (9) can be found by introducing a Green's function. The final 
result for (on the surface) is then 

4>°(x,u) = F[(e h -v h + e z K)e ], (10) 

where Vh is the horizontal gradient. The zero-order surface velocity potential is given by 
e°(x, z = 0, u) = 2vr J J G(x - x', u)S°(x', cu)d 2 x', (11) 

with 

G ^=(k°JJ^ir/ k - (12) 

where k = \\k\\ is the horizontal wavenumber. In summary, equations (9)-(12) give the 
observed zero-order filtered wave field at the surface for f modes propagating through the 
background model with no flows present. 



3. Travel-time shift due to a uniform flow 

Here we study the effect of a horizontal, spatially-uniform steady flow, 

u = u x x + u y y = const. (13) 

From equations (2)-(5) we see that the effect of the flow is to transform u into uo — k ■ u. 
The Doppler shift k ■ u is the result of a simple Galilean transformation x — > x — ut. This 
transformation applies to all quantities in the problem except the filter F, which is applied 
independently of the flow. The filtered data cube can thus be expressed in terms of the 
Doppler velocity in the absence of a flow, 0°, according to 

tj}{k, u) = F(k, uj)<j)°(k, u-k-u). (14) 
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The fundamental computation in time-distance helioseismology is the cross- covariance 
of the oscillation signal, measured for two different points x± and x 2 on the solar surface, 

C(cci, x 2 , u) = — ip*(x 1 , u)if){x 2 , u), (15) 

where the star denotes complex conjugation and T is the observation duration. Because of 
horizontal translation invariance (uniform flow), the expectation value of the cross- covariance 
is the inverse Fourier transform of the power spectrum (e.g. Gizon & Birch 2004): 

C(x 1 ,x 2 ,u) = JJ P(k,u)e ik A d 2 k, (16) 

where A = x 2 — x\ and 

P(k,io) = EMk,io)\ 2 ] (17) 

is the expectation value of the power spectrum of oscillations. The computation of P requires 
the knowledge of the source covariance function m s (u); here we choose the same function 
as in Gizon & Birch (2002). A general configuration of the observation points and notation 
used throughout the paper is depicted in Figure 1. 

We now consider two strategies for the computation of the cross-covariance C in the 
case of a uniform flow. The first is to compute equations (16) and (17) numerically, to obtain 
what we call the exact solution. The second is to compute successive approximations to C 
using a Taylor expansion in u of the power spectrum P. For a constant flow, we have 

P{k,u) = P Q {k,u)-k-u\F\ 2 d u [\<f\ 2 ] 

+ \{k.uf\F\ 2 di[\<t>r] 

- i(fc.«) s |F| 2 ^[I^T]+---, (18) 

where P°(k,uj) = E[\ip°\ 2 ] is the expected power spectrum when u — 0. We may truncate 
this expansion after each successive term to study the zero-, first-, second-, and third-order 
cross-covariance (terms in u°, u 1 , u 2 and u 3 respectively). 

Figure 2 shows the cross-covariance function at each level of approximation of the Taylor 
expansion for a uniform flow. Convergence to the exact covariance is observed as higher- 
order terms are retained. Notice that the phase shift is largely accomplished by the first-order 
term, while the amplitude correction is largely accomplished by the second-order term. 

In this paper, we are mainly interested in studying travel-time differences, which are 
obtained by finding the time t + (x 1 ,x 2 ) it takes for waves travelling from x x to x 2 , and 
subtracting the time T-(x 2 , Xi) for those waves going from x 2 to X\. In other words, 

T diff (cci, x 2 ) = t + (x 1 , x 2 ) - r_(£c 2 , x-s). (19) 
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This quantity gives information about the underlying flows. Following Gizon & Birch (2004), 
the travel-time difference can be measured from the cross-covariance function according to 

2 Im J°° uo C°(A, uo) AC(x 1: x 2 , uo) duo 

x 2 ) = j^Wi^W^' ' (20) 

where C°(A,u;) = C°(x 1 ,x 2 ,uo) is the covariance function in the absence of a flow and 
AC = C — C° is the change in the cross-covariance due to the flow. The operator Im takes 
the imaginary part of the expression. 

Using the zero-, first-, and third-order cross-covariances, as well as the exact cross- 
covariance, we now study travel-time differences using equation (20). This will provide us 
with a means to test the first-order Born approximation which is developed in § 4.1, as well 
as to quantify higher-order terms. Note that the second-order term does not affect travel- 
time differences because it only introduces perturbations to the cross-covariance that are 
symmetric in the time lag. 

We first study the travel-time differences for various separation distances A and flow 
speeds u. Figure 3 shows the exact and approximate scaled travel-time differences r^s/it 
versus A computed at fixed u = mA for u = 200 ms 4 and u = 400 ms _1 (A = A/ A). We 
find that the first-order approximation is within 10% of the exact value for u = 200 ms 4 
and A < 20 Mm, while the error is only about 1% when terms up to u 3 are kept. For a flow 
that is twice as large, u = 400 ms" 1 , non-linear effects are more important. The functional 
dependence at large A is partly due to mode damping, as will be shown later in § 4.7. 

Next we compare this result with observed travel times. We used 4 hours of quiet Sun 
MDI Dopplergrams tracked at a constant rate of 678 ms" 1 with respect to the local surface 
rotation rate in order to mimic the effect of a constant flow. Travel times are then measured 
using equation (20). Figure 4 shows the result. We find a good agreement between the exact 
model and the observations. The agreement is not, however, perfect at short distances; this 
is likely to be due to a mismatch between the model and the solar f-mode power spectra. 

Figure 5 shows the variations of r diff as a function of u for fixed distances A = 5 Mm 
and A = 10 Mm. The first-order approximation is within 10 % of the exact solution for flows 
with amplitudes less than about 400 ms -1 , although for larger A the agreement becomes 
worse. However, these results demonstrate that the first-order approximation (and most 
likely the kernels given in § 4) is appropriate for the study of supergranular flows. 

Let us turn to the question of whether advection of the waves by the flow is the only 
important effect. To first order in u, the change in the cross-covariance in the case £ = z 
has three distinct components: 

AC ~ -(2tt) 6 [[ d 2 ke ik -*(k-u)\F\ 2 
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x [m 8 |K| 2 9 w (|G| 2 ) + (oLm a )|«| 2 |G| 2 

+ m s |G| 2 ^(|/€| 2 )] . (21) 

The three terms in brackets in this equation correspond to the effect of the flow on the waves, 
on the sources, and on the damping. We can thus determine the relative contribution that 
each of these physical effects separately has on the first-order travel-time differences. This 
is shown in Figure 6. The advection of the waves contributes by far the most to the travel- 
time difference, making up about 98% of the total value. We conclude that, for practical 
applications, we can ignore the effects of the flow on the damping and source functions. 



4. Spatially-varying horizontal flow 

Now we relax the assumption of a constant flow and consider a spatially- varying, steady, 
horizontal, small-amplitude flow u(x). The general study of the interaction of surface waves 
with flows involves a number of complications, and we do not attempt to solve the problem in 
its full generality; instead we require the flow vorticity be zero. Due to the small amplitude 
of the flow, the problem may be solved in the first Born approximation, which we carry out 
in § 4.1 and will lead to the sensitivity kernels of § 4.2. 



4.1. First-order Born approximation 

Here we solve the linearized set of hydrodynamic equations introduced in § 2 by applying 
the first Born approximation. We expand each wave quantity q into a zero-order term, q°, and 
a perturbation term due to the flow, 5q. We then insert the expanded wave quantities into 
equations (2)- (5) and retain only terms which are of first order. In the bulk, the perturbed 
velocity potential <56 satisfies the three-dimensional Laplace's equation, 

V 2 56 = 0. (22) 

Eliminating the pressure and the surface elevation, the surface boundary conditions at z = 
reduce to 

(d z - k)5Q = iu ■ V [(d w /c)e° + d^S ] = 5S. (23) 

In doing so, we have used the approximations 5U ~ iu ■ Vc^n and 5TQ° ps 1(8^)11 ■ V@°, 
which describe the perturbations to the wave sources and the wave damping respectively. 
We note that SO satisfies the same equations as 0°, except for a different source function 
SS on the right-hand side of the surface boundary condition, equation (23). This is the first 
Born approximation, an equivalent-source description of wave scattering. 
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Using the same Green's function as in § 3, the first-order surface velocity potential is 

SG(x, z = 0,u)=2n JJ G(x - x' , cu)5S(x', u)d 2 x', (24) 

from which we deduce the perturbation to the filtered data, 

6ip(x, u) = F [(4 • V h + £ x k)6G + ii z {d w n)u • V0°] . (25) 

4.2. Linear sensitivity kernels for flows 

The first-order perturbation to the cross-covariance is obtained by approximating AC 
in equation (20) according to 

AC(jei, x 2 , w) ~ E [5^*(x 1 ,u)^°(x 2 ,uj) 

+ V°*(aJi,w)^(aj2,w)], (26) 

where 5ip is the first-order perturbation to ip. Using the expressions we have found for ip 
and 5ip, equation (26) can be put into the form 

AC(xi, x 2 , oS) ~ JJ C(x 1 , x 2 , co] x) ■ u(x) d 2 x, (27) 

where the function C, explicitly written in Appendix A, gives the linear sensitivity of the 
cross-covariance to a spatially-varying flow. Inserting equation (27) into equation (20) gives 

rdifffci, x 2 ) = JJ K diS ( Xl , x 2 ; x) ■ u(x) d 2 x, (28) 

where the function K dlfr is the two-dimensional sensitivity kernel, defined by 

2 Im r o °°^C (A,^)C(a; 1 ,a; 2 ,^;a;) d^ 

J °°uj' 2 \C°(A,uj')\ 2 duj> ' 1 ' 

Figure 7 shows example kernels for the linear sensitivity of the travel-time difference to 
small- amplitude flows, for a separation of 10 Mm between the observation points and for a 
vertical line of sight £ — z. A short discussion of the numerical computation of kernels is 
given in Appendix B. The kernel in Figure 7a, K dlff , gives the sensitivity of the travel times 
to u x and the kernel in Figure 7b, K dlff , gives the sensitivity to u y . For reasons of symmetry 
(the Ky lS kernel is antisymmetric with respect to the lines x = and y — 0) the total integral 
of Ky lS is zero. 

We observe elliptical and hyperbolic features in the kernels, similar to what was seen in 
the kernels for the damping rate and the source strength perturbations derived by Gizon & 
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Birch (2002). The elliptical features have been called Fresnel zones in geophysics (e.g., Tong 
et al. 1998). The first Fresnel zone bounds the large area between the observation points in 
Figure 7 where the kernel is nee; ative. That Kf s is negative in this region is intuitive: local 
flows in the +x direction should speed up waves traveling in this direction while slowing down 
waves going the opposite way, resulting in negative travel-time differences. The hyperbolae 
and ellipses, due to scattering of waves generated by distant excitation events, are discussed 
futher in § 4.6. 

It is important and instructive to check the Born kernels for consistency with the first- 
order approximation of t&s/u calculated under the assumption of a constant flow in § 3. It 
also enables us to verify the numerical computation of the kernels. If we again consider a 
uniform flow, u = ux, then from equation (28) it is possible to spatially integrate the kernels 
to obtain 

K* ff ( A; x) d 2 x. (30) 

We have computed kernels for several values of A and obtained the corresponding t&s/u. 
The results are shown as large dots in Figure 3. The total integral of is within 0.1% of 
the value expected from the first-order approximation, giving us confidence in the numerics. 




4.3. Effect of the line-of-sight projection 

The kernels in Figure 7 were computed with the line-of-sight vector £ = z, i.e., at solar 
disk center. We would like to obtain kernels at any general point on the solar disk. The 
variation in the line-of-sight vector across the solar disk is incorporated in our model and 
appears explicitly in the expression for the kernels (see Appendix A). Foreshortening in the 
center-to-limb direction (the point-spread function effectively depends on the direction of k) 
is another important systematic effect, which has been ignored for the sake of simplicity. 

A region on the Sun around the solar longitude measured from disk center, $, and 
latitude, L, has a local line-of-sight vector approximately given by 

£ = — sin $ x — sin L cos <3> y + cos L cos $ z. (31) 

Two examples of kernels computed away from disk center are given in Figure 8 for coordinates 
$ = 7r/3 and L = (on the equator towards the western limb). The kernel in Figure 8a 
is K^ lff , whereas the kernel in Figure 8b is K^ lff , whose observation points lie on the same 
meridian (A = y). If the line-of-sight effect did not matter, then the two kernels in Figure 8 
and the kernel of Figure 7a would all be identical after rotation of Figure 8b. This is not the 
case, however. 
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In particular, although these three kernels have the same total integral (see below), the 
kernel of Figure 8b looks very different. This is because, in the case of Figure 8b, waves 
that scatter along the line that connects the two observation points have a small component 
along the line of sight and have thus a small contribution to the cross-covariance function. 
In general, the flow sensitivity is reduced when the scattered wave at the observation points 
is travelling perpendicular to the center-to-limb direction. For example, there is about 50 % 
less sensitivity along the ray path connecting the two observation points for the kernel of 
Figure 8b compared to the kernel of Figure 7a at disk center. 

We have computed the total integral of the kernel as a function of position on the solar 
disk. We studied two cases with L = and varying <3> in the range from — it/2 to 7r/2: one 
in which the observation points lie on the equator and the other in which the observation 
points lie on the same meridian. The resulting travel-time differences show very little, if 
any, variation across the disk. This means that large scale flows are practically insensitive to 
line-of-sight effects. However, judging by the fine structure of the kernels, the line-of-sight 
effect is important for small-scale flows. 

4.4. Source advection is unimportant 

In the uniform flow model it was possible to separate the effects of advection of the 
waves, of the wave sources, and of the wave damping to determine the contribution of each 
to the travel-time difference (see eq. [21] and Fig. 6). We concluded that for uniform flows the 
wave advection is by far the dominant contribution to the travel time shifts. Is this also true 
for a spatially- varying flow? Using equation (A3), it is possible to identify the contribution 
to the kernel due to source advection. At A = 10 Mm we find that the source advection 
contribution is about three orders of magnitude less than the contribution from the advection 
of waves by the flow. This result is useful for future kernel calculations. Neglecting source 
advection is an excellent approximation. This approximation has been used in the past by, 
e.g., Woodard (2002), Birch & Gizon (2007), and Birch et al. (2007), although not explicitly 
justified. 

4.5. Kernels depend on central frequency and bandwidth 

Birch et al. (2004) showed how phase-speed filtering affects the properties of p-mode 
sensitivity kernels for sound speed perturbations. Here we study how f-mode kernels change 
as the spatial frequency content of the input f-mode power spectrum is varied. To isolate 
particular f modes, we apply a multiplicative filter G{k) to the wavefield, in addition to 
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the filter F that contains the instrumental MTF and a mode- mass correction (§ 3). This 
additional multiplicative filter is a simple Gaussian function, 

(32) 

with central wavenumber k and standard deviation a. Note that f modes with different 
wavenumbers probe different depths (e.g. Duvall & Gizon 2000). 

Figure 9 shows examples of sensitivity kernels, using the filter G(k) for k R = 487, 
900, and 1322 and aR = 348 (R = 696 Mm is the radius of the Sun). Plotted alongside 
each kernel is the corresponding filtered power spectrum. The particular part of the f-mode 
ridge that is isolated by the additional filtering is evident from these power spectra. As the 
bandwidth is reduced by filtering, the kernels show more fine structures and are more spread 
out in space (see Figure 7 for a comparison with a kernel computed with the full f-mode 
power). As ko is varied from smaller to larger values, more and more fine structure becomes 
apparent in the kernels. We conclude that the sensitivity kernels must be tuned to account 
for the correct power spectrum before any high spatial resolution inversion is attempted, 
especially at scales comparable to the wavelength. 

4.6. Far-field approximation to the kernels 

We now consider an approximation to the sensitivity kernels given by equation (29). 
This is carried out in the far field, i.e., for points x away from the observation points a?i and 
X2- The solution will enable us to speed up the numerical computation of the kernels, as well 
as to expose some of the main features of the kernels that have already been discussed. The 
approximation allows for the analytical solution of the integrals given in equations (A6)-(A9). 
The details are presented in Appendix C, where the final expression for the approximate 
kernel is equation (C8). 

Figure 10 is an example sensitivity kernel computed in the far-field approximation. It 
resembles quite well the 'exact' kernels of § 4.2 (see Fig. 7a). The agreement between the 
exact and far-field kernel is everywhere within 10 % except within a 1 — 2 Mm radius about 
the two observation points. In fact, the far-field kernel is not defined at the observation 
points. Figure 11 shows one-dimensional cuts along the x = line of both the exact and 
approximate K^ lfr . 

The far-field approximation kernel, equation (C8), clearly brings out two important 
features. The first is the origin of the ellipses and hyperbolae mentioned in § 4.2. Both 
the ellipses (Ai + A 2 = X = const) and hyperbolae (Ai — A 2 = A = const) are clearly 
present in the numerator of expression (C8). The second is the presence of damping through 



G(k) = exp 



-(k - k y 

2a 2 
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the K\ terms (k{ = Im(/t) = ju/g). The damping causes the kernel to decay by a factor 
~ exp[— uj~f(Yl + A)/g] as one moves away from the ray path, where 7 is the damping rate. 



4.7. Comments on the ray approximation 

Ray theory is a high-frequency approximation. According to Kosovichev & Duvall 
(1997), in the ray approximation the sensitivity of travel times to flows is confined to the 
ray path which joins the two observation points. In the case of f modes, the ray path is a 
straight line between the two observation points. For simplicity we study a configuration 
where the two observation points lie on the x axis, cci = (0, 0) and x 2 = (A, 0). The one-way 
travel-time perturbation becomes 

1 ' A 



S T ™y = — / 5k x (x)dx, (33) 
^0 Jo 

where 5k x is the local change in x component of the wavevector and u = 3 mHz is the 
frequency for which the /c-averaged power spectrum is maximum. For the case of f modes 
and a flow u x (x), the travel-time difference t^s can be approximated as 

tZ = ~ A 4- f«*te. (34) 
9 Jo 

We would like to check whether the two-dimensional kernels derived in this paper reduce 
to the ray kernels when we compute the cross-path integral, i.e., when flows vary slowly in 
the y coordinate. In order to facilitate an analytic expression for the travel times, we begin 
from the far-field approximation of K^ lff , equation (C8). Following Tong et al. (1998), the 
integral of K^ lff over y can be estimated by the method of stationary phase (sp). We obtain 
the following 

[lV f u x dx, (35) 
Jo 



T, 



sp 



with the definition 



f-pye-^du, 



1/2 



"(A) - 1^2 f o °°p2 e -2 Ki A cos 2 (/s . rA _ ^4)^1 ' (;J)()) 

where P^(cu) is an average of the power spectrum over wavenumbers, and is explicitly defined 
in Appendix C. The functions k t (cu) and Ki(u) are the real and imaginary parts of the 
complex wavenumber at resonance. For A > 5 Mm, we find that ZJ ~ 3.3 mHz is a weak 
function of A. At shorter distances however, uJ(A) varies fast with A. Only in the case when 
the power distribution P^(u) peaks around uq, at large distances A, and in the absence of 
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damping, do we have uJ = u and thus recover the ray approximation given by equation (34). 
In practice, damping cannot be ignored, and the naive ray kernel should be corrected by 
replacing u by uJ(A). 

Figure 12 shows the travel-time differences for a uniform flow from equation (35) for two 
cases: one with damping present and one with damping set to zero. The difference between 
the two cases is quite large. Damping reduces the overall sensitivity to the flow (thick solid 
line). The nonlinearity of the travel times with distance is a direct effect of the damping, 
and we see that turning off the damping restores linearity (dot-dashed line). For comparison 
we also plot the exact travel-time differences for a uniform flow of 200 ms" 1 , studied in 
§ 3. Remarkably, for A < 20 Mm, the exact solution and the solution from equation (35) 
are in good agreement, even though there have been two approximations made to arrive 
at equation (35). The near-field wiggle at ~ 2.5 Mm is captured by the improved ray 
approximation. 



4.8. Sensitivity kernels for mean travel times 

Let us consider the mean travel times for point-to-point measurements defined as 

Tnwanfcl, X 2 ) = ^ [r+(x 1 , X 2 ) + T_(x 2 , SEi)] , (37) 

where r± has been previously defined. It is sometimes assumed that mean travel times are 
not sensitive to flows. This is not strictly true, as we now discuss. 

Following Gizon & Birch (2002) and § 4.2, we may derive linear sensitivity kernels for 
mean travel times which satisfy 

Tmeanfcl, X 2 ) = JJ K™^, X 2 , x) ■ u(x) d 2 X, (38) 

where 

_ -2Re / °° Hilb [uC°(A, hj] C(x lt x 2 , U ; x) du 

r o °°o;' 2 |CO(A,a;0| 2 cla;' ' 1 } 

Here, Re takes the real part of the expression in the numerator and Hilb[. . .] denotes the 
Hilbert transform. 

To show how spatially-varying flows are related to mean travel times r mea n, we compute 
the corresponding sensitivity kernels K mean and plot them in Figure 13. First, we point 
out that the spatial integrals of K™ ean and K™ ean are both zero, indicating that a constant 
flow has no first-order effect on mean travel times. There are, however, some structures of 
significant magnitude in these kernels on small length scales (as large as for the travel-time 
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difference kernel K dlff ). Thus, flows that vary on a scale comparable to the wavelength may- 
leave their signature in the mean travel times. 



Birch & Gizon (2007) computed three-dimensional kernels for the effect of local flows on 
travel-time differences, taking into account solar stratification but neglecting the advection 
of sources and damping. As we have shown here, the effects of source and damping advection 
are, for practical purposes, negligible compared to the direct advection of the waves by the 
flow. Hence, it is meaningful to compare the two calculations. 

The 3D kernels, denoted here as K 3d , give the linear sensitivity of travel-time differences 
to general three-dimensional steady flows, u(x,z), and satisfy 



In order to compare the three-dimensional f-mode kernels, K , with the two-dimensional 
kernels computed in the previous sections, we define the depth-integrated kernel 



where the integral is over all z. The kernels K BG give the sensitivity of travel-time differences 
to flows which are independent of depth. 

Figure 14 compares an example depth-integrated kernel, K BG , with the corresponding 
two-dimensional kernel from this paper. The total integrals of the two kernels agree to within 
a few percent. Notice, however, that the details of the fine structure are somewhat different. 
These differences are likely the result of the somewhat different assumptions regarding the 
source covariance, the damping model, and the instrumental point-spread function. 

The two bottom panels of Figure 14 show one-dimensional cuts through the two kernels. 
The sizes of the first Fresnel zones are equal; this shows that the dominant wavelengths in 
the two calculations are similar. The kernel K BG shows somewhat less ringing than the 
kernel computed in this paper, which could be due to a slightly different frequency content. 
The differences in the fine structure of these two kernels demonstrate the importance of 
accurately modeling the true bandwidth of solar f modes (see § 4.5). 



5. Comparison with simplified three-dimensional kernels 




(40) 




(41) 
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6. Discussion 



We have shown how to compute f-mode travel-time shifts due to uniform flows. The 
result depends on travel distance in a complicated way for several reasons: near-field effects 
at short travel distances and damping at large distances. In addition, we find that the 
linearization of the travel times with the flow amplitude is valid to within 10% for flows 
with amplitudes less than about 250 ms" 1 and travel distances less than about 25 Mm. 
For larger flows or larger distances, non-linear effects become important, as the travel times 
do not scale like the flow. This non linearity may complicate the interpretation of time- 
distance measurements because large-amplitude flows are observed on the Sun, including, 
e.g., rotation (hence the need to at least remove the main effect of rotation before applying 
time-distance helioseismology) . 

We computed the two-dimensional sensitivity of f-mode travel times to spatially- varying 
steady flows, using the first Born approximation (small flow regime). We found that the main 
physical mechanism responsible for travel-time shifts is, by far, the advection of the waves 
by the flow. The effects of the advection of wave sources (granulation) and wave damping 
by the flow are completely negligible. This will substantially simplify future computations 
of kernels, and justifies the omission of these terms in previous calculations. The kernels are 
also highly sensitive to the frequency content (central frequency, bandwidth) of the input 
f-mode power. Correctly tuning the model power spectrum to the observed one is therefore 
crucial for the use of these kernels in any inversion. The point-to-point kernels computed here 
could easily be averaged to obtain point-to-quadrant or point-to-annulus kernels, a standard 
averaging technique used in time-distance helioseismology. 

We showed that, at small spatial scales, the sensitivity kernels for Doppler velocity mea- 
surements depend very significantly on position on the solar disk and on the angle between 
the observation points and the center-to-limb direction. This seriously complicates the inter- 
pretation of travel-time measurements as the inversion problem is not a deconvolution. We 
note that close to the limb, kernels are nothing like ray theory would suggest. We caution 
that we have not considered the effects of foreshortening in this paper (see Jackiewicz et al. 
2006, for a preliminary discussion). 

The two-dimensional f-mode kernels are fast to compute (especially in the far field) and 
compare well with depth-integrated three-dimensional kernels. They may prove quite useful 
to infer near-surface horizontal flows. 
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A. C function 



We would like to derive an expression for the function C in equation (27): 

AC(x 1 ,x 2 ,u) = J J C(x 1 , x 2 , ou; x) ■ u(x) d 2 x. (Al) 

First, we use the zero-order and scattered wave fields in equations (10) and (25) and plug 
them into the first-order approximation to the cross-correlation, given by 

AC(x 1 ,x 2 ,ou) « E [5iP*{x 1 ,uj)iIj (x 2 ,uj) + iP°*(x 1 ,uj)5iIj{x 2 ,uj)] , (A2) 

where the flow u is contained in the 5ip terms. After some manipulation one can show upon 
comparison of equation (Al) and (A2) that 

C( Xl , x 2 , u; x) = im s [(cU)!T(A 2 , u) V X I(A 1: u) - (d^*)K(A u u;)V x I*(A 2: lu)] 

+ i(oLm s ) [K*(A 2 , uj)V x K(A 17 uj) - K(A 1 ,u)V x K*(A 2 , oo)} , (A3) 

where we use the notation Ai = x — xi, A 2 = x — x 2 , and A = x 2 — x\ (see Fig. 1 for the 
geometry schematic). The function m s is equivalent to the source function in equation (9) 
for spatially uncorrelated sources, and is related to the m° in equation (51) of Gizon & Birch 
(2002), where m s = uj 2 m /g 2 . In the above equation, we defined the following functions: 



where 

Jo(x,u) 
Jo(x,u>) 

Jl(x,Ld) 

Ji(x,u) 

The second line in equation (A3) is the Doppler shift, or advection, of the source function 
due to the flow. This term is small compared to the wave advection term in the first line, 
as discussed in § 4.4. We note that a complicated term has been left out of equation (A3) 
that contributes only at the observation points X\ and x 2 . We have calculated it, yet it is 
negligeable for practical applications and so we have ignored it. 



/ = t z Jo + £ h -Ji, (A4) 

K = £ z J + i h -Ji, (A5) 

= (27r)Vjy e ik - x \G(k,uj)\ 2 F*(k,uj)d 2 k, (A6) 

= (27r)Vjy e ik - x G*(k,uj)F*(k,uj)d 2 k, (A7) 

= -(2vr) 5 i J J e ik - x k\G(k,Lu)\ 2 F*(k,iu)d 2 k, (A8) 

= -(2vr) 2 i [[ e ik - x kG*(k,u)F*(k,uj)d 2 k. (A9) 
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B. Numerical computation of kernels 



We briefly point out some aspects of the numerical computation of the kernels. The 
most important step which determines the length of the computation is setting up the initial 
x — y grid on which the kernels will be calculated. The kernels of Figure 7, which lie on a 
horizontal spatial grid of resolution 0.2 Mm in both x and y directions (145 x 145 points), are 
computed in ps 30 minutes on a single processor machine, for a code written in C++. However, 
this value of the resolution turns out to be unnecessary. One may use a much coarser grid, 
and then perform a two-dimensional Fourier interpolation to recover the kernels onto a finer 
grid, if necessary. We have tested this and found that the coarsest grid we may use for these 
kernels is roughly of a resolution of ~ 1 Mm. The Fourier interpolation method is due to 
the Nyquist-Shannon sampling theorem. This theorem tells us that the minimum resolution 
necessary to recover the full signal is da; = A/4, where A is the dominant wavelength in the 
problem, which, for f modes at power maximum, is about 5 Mm. Implimentation of this 
speeds up the computation quite effectively The spacing for the k and u grids used for the 
computation are approximately duo = 7 x 10~ 5 rad s _1 and dk = 0.01 Mm -1 , over the ranges 
oj = [1.5,5] mHz and k = [0,2.5] Mm' 1 . 



C. Far-field approximation to the kernels 



We would like to find an analytic approximation to the expression for the sensitivity 
kernels, starting from equation (29): 

K-diff _ ~ 2Re JiT mC°(A, u) C{x x , x 2 ,uj; x)duj 

r o °°a;' 2 |C o (A,a;0|W ' 1 j 

We consider for simplicity sensitivity kernels where £ — z, and we ignore the source advection 
term in C (line 2 of eq. [A3]) , which we have already shown to be quite small. 

Using the definition of the power spectrum (eq. [17]), the zero-order cross-correlation 
(eq. [16]) can be written as 

l^^)| 2 „i fe .A j2 



C u (A,cu) =m 8 {u)\K\ 



\k — k\ 



d z k, 



(C2) 



where m s (u) is the source function described in Appendix A, F is the filter function, and n 
is the complex wavenumber. Similarly, C of equation (A3) can be expressed as 

12 rr Ud^)F(k,u)F*(k',cu) ^ ik . A2+ik , Ai 
k\ 2 k — K 



im s |ft| 2 ff i/ 
7 ~{2^JJ |F= 



(d^)F*(k,u)F(k',u) t 
k — k* 



,ifc-Ai-ifc'-A 2 



d 2 kd 2 k'. (C3) 
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The integrals over wavenumber in equation (C2) and (C3) can be performed analytically 
when the argument of the exponential functions is large, i.e., in the far-field limit. For any 
smooth function f(k) and horizontal vector x, we have, in the limit of large x — \\x\\ : 

t^— tfk » fi^^ f{n x x)e^ x -^e-^ (C4) 



a: 

JJ-oo \k-K\ z 2lKi J J k-K ZlKiJJ k-K* 



» (27r)3/2 ^l Ki:C [f(K T x)e [ ^ x ~^ + f(-K T x)e^ x -^] , (C5) 

where n r = Re(/c) = w 2 /g and Ki = Im(/t) — ^fuj/g (see § 2). 

Using these integral approximations we can now evaluate the kernel in equation (CI). 
For simplicity, let the filter function be real and independent of k, and to simplify the 
notation let F(k, u) — > F. Integrating equation (C2) with the help of (C5) yields 

C°(A,uj) ~ Pv> !_ I A cos(/t r A - tt/4), (C6) 

where P^{uj) = (2Tc) 3 ^ 2 m s F 2 (K r , uj)k 3 /ki is approximately the power spectrum integrated 
over wavenumber. Note that this term is real. Similarly, from equation (C3) one can show 
that 

Re[i^C] = s ] [cos(/t r S - tt/2) - cos(« r A)] ( A 1 - A 2 ) , (C7) 

/CiVAiA 2 V / 

where £ = Ai + A2 and A = A2 — Ai, A, denotes the unit vector in the direction of 
Aj, and we have used \k\ 2 ~ k 2 since the imaginary part is small (damping). Finally, 
plugging equations (C6) and (C7) into equation (CI) gives the far-field approximation to 
the sensitivity kernel: 

-2 (Ai - A 2 ) p^-^+A) [ cos ( KrS _ n / 2 ) - cos k t A] cos(/t r A - 7r/4)dcu 

R1 {Xl,X2 ' ,X) *~JW^/X io°° P 2 e~ 2K[A cos 2 (/t r A — 7r/ 4)d^j ' 

(C8) 

As noted in § 4.6, the kernel is not defined at the observation points (A x , A 2 = 0). A plot 
showing this kernel is given in Figure 10. Expression (C8) reduces to the result of ray theory 
in the infinite frequency limit, as demonstrated in § 4.7. 
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Fig. 1. — Schematic of the geometry of the plane-parallel patch of the Sun we consider. 
Travel times are measured using the wave signal at the two photospheric observation points 
Xx and x 2 - The distance between these points is A = ||A||. The vectors given by A 1 and 
A 2 appear in the expression for the sensitivity kernels (see Appendix A). 
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Fig. 2. — Temporal cross-covariance functions C(xi,X2,t) in the uniform flow model for 
u = 400 ms 4 A and A = 10 Mm. The thin solid lines are the exact cross- covariances 
(repeated in each panel), and the thick dashed lines correspond to the approximate cross- 
covariances, truncated at the order indicated in the legend of each panel. 
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A (Mm) 

Fig. 3. — Scaled travel-time differences in the uniform flow model as a function of separation 
distance A for two different flow speeds, where u = ua. The thick solid line is the linear 
approximation (which is independent of the flow when plotted as r/u). The upper and 
lower dashed lines correspond to the exact solutions for u = 400 ms" 1 and u = 200 ms" 1 
respectively. The thin solid lines are the third-order approximation for each u. The big dots 
denote the travel-time differences derived from integration of the kernel Kf s discussed in 
detail in § 4.2. 
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Fig. 4. — Comparison of travel-time differences from data and from the exact solution for 
a uniform flow. The circles are the travel-time differences measured from quiet sun MDI 
Dopplergrams tracked at the constant rate of 678 ms _1 . The solid line denotes the travel- 
time differences computed from the exact model of § 3 for a constant flow of 678 ms 4 . The 
one standard deviation error bars of the data are also shown. 
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Fig. 5. — Travel-time differences as a function of uniform flow speed u for two different 
distances A. For each set of curves, the dashed line is the exact solution, the thick solid 
line is the first-order approximation, and the thin solid line is the third-order approximation. 
The top set (top three curves) of travel-time differences is for A = 5 Mm, while the bottom 
set (bottom three curves) is for A = 10 Mm. 
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Fig. 6. — First-order travel-time difference contributions from the three separate advection 
terms of equation (21) for a uniform flow. The y coordinate scale is the log of the absolute 
value of the travel-time difference per second. In this case, A = 10 Mm and u = ma. 
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Fig. 7. — F-mode travel-time difference sensitivity kernel for flows as a function of horizontal 
position x = (x,y). Panel (a) shows K^ lff and panel (b) shows K^ lff , with A = 10 Mm x 
in both cases. The crosses denote the two observation points at x-i = (—5, 0) Mm and 
X2 = (5,0) Mm. The units of the kernel are sMm~ 2 (km/s) -1 . 
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Fig. 8. — F-mode travel-time sensitivity kernels computed away from disk center, at $ = 7r/3, 
L = (on the equator towards western limb). The definition of the angles is given by 
equation (31). The two observation points are denoted by the crosses, (a) K^ lff for a?i = 
(-5, 0) Mm, x 2 = (5, 0) Mm. (b) Kf s for the case when x 1 = (0, -5) Mm, x 2 = (0, 5) Mm. 
The direction towards the disk center is annotated. The units are sMm -2 (km/s) -1 . 
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Fig. 9. — Power spectra (a), (c), and (e) and sensitivity kernels Kj ff (&), (d), and (/) 
computed using the additional filter G(k) defined in equation (32). The filter used for the 
top row has a central wavenumber koR = 487, for the middle row k$R = 900, and for the 
bottom row k^R = 1322. In each plot aR = 348, and R = 696 Mm. For each kernel, 
A = 10 Mm, and the units are sMm -2 (km/s) -1 . The dispersion relation for f modes, 
uj 2 = gk, is given by the dashed lines in the power spectra. The gray scale of the power 
spectra is linear and truncated for emphasis. 
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Fig. 10. — Far-field kernel K^ lfr computed from equation (C8) and discussed in § 4.6. Note 
that this kernel is not defined at the two observation points, given by the crosses. The spatial 
integral of this kernel is within 5 % of the spatial integral of the exact kernel (see Fig. 7a). 
The units are sMm -2 (km/s) -1 . It is computed in just several seconds. 
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Fig. 11. — Cuts through the exact and far-field kernel along x — 0. The heavy solid line 
is the cut through the kernel of Figure 7a, and the thin line is a cut through the far-field 
approximate kernel from Figure 10. 
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Fig. 12. — Scaled travel-time differences in the ray approximation, the stationary-phase 
approximation, and for the exact model. For each case, we consider a constant flow in 
the direction of x x to x 2 - The exact solution {dashed line) is for a flow u = 200 ms" 1 
(see § 3). The thin solid line corresponds to travel-time differences computed using 
equation (34). The thick solid line is computed from equation (35) and denotes t^ s . Finally, 
the dot-dashed line is the travel-time difference from equation (35) without any damping. 
The central angular frequency ujq = 3 mHz. 
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Fig. 13. — Sensitivity kernels for mean travel times r mean discussed in § 4.8, where (a) is 
K~ and (b) is K~ and A = 10 Mm. The units are sMm~ 2 (km/s)" 1 . 
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Fig. 14. — Comparison of the x components of K BG (the depth-integrated kernel from Birch 
& Gizon 2007) and the corresponding two-dimensional f-mode kernel computed in this paper, 
K difr , for the case A = 10 Mm. Panel (a) shows K BG . Panel (b) shows K dlff . The horizontal 
integrals of the two kernels agree to within a few percent. Panels (c) and (d) show one- 
dimensional cuts through the K BG kernel (thick lines) and the two-dimensional kernel K dlff 
(thin lines) along the y = and x = lines respectively. In all panels the units of the kernels 
are sMm~ 2 (km/s) -1 . 
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